function A = uf(A)
% Computes the following
% uf(A) = A/sqrt(A'*A);
% 
% This extracts the orothogonal factor of the polar decomposition of A
% 

[L, ~, R] = svd(A,0);
 A = L*R';
